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The superconducting pairing of electrons in doped graphene due to in-plane and out- 
of-plane phonons is considered. It is shown that the structure of the order parameter 
in the valley space substantially affects conditions of the pairing. Electron-hole pair- 
ing in graphene bilayer in the strong coupling regime is also considered. Taking into 
^^ I account retardation of the screened Coulomb pairing potential shows a significant 

CN ■ competition between the electron-hole direct attraction and their repulsion due to 

virtual plasmons and single-particle excitations. 
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c/3 , 1. Introduction 

Electrons in graphene, a two-dimensional form of carbon, can be described by a 
two-dimensional Dirac-type equation for massless particles near the Fermi level (see 
C^ ' Castro Neto et al. (2009) and references therein). Thereby, graphene offers a unique 

possibility to study effectively ultrarelativistic charged particles in condensed mat- 
, ter phenomena (Katsnelson et al. 2006, Katsnelson & Novoselov 2007) and particu- 

^ . larly in collective phenomena (Lozovik et al. 2008, Berman et al. 2008a,&). Absence 

O ' of mass for electrons make it possible to achieve new regimes of quantum many- 

O . particle systems behavior in graphene. Therefore it is interesting to search for vari- 

ous superconducting and superfluid phases in graphene and graphene-based struc- 
tures, with their applications for dissipationless information transfer in nanoscale 
^ ' devices. In the present paper, we consider the Bardeen-Cooper-Schrieffer-like (BCS- 

t^^ , like) (Bardeen et al. 1957) phonon-mediated pairing of electrons in graphene, and 

the Coulomb pairing of electrons and holes in graphene bilayer, taking into account 
the unusual electron dynamics. 

Several possibilities for electron pairing phenomena in graphene were proposed. 
One possibility is the pairing of spontaneously created electrons in the conduc- 
^^ ' tion band and holes in the valence band, leading to the excitonic insulator state 

(Khveshchenko 2001). The results of numerical simulations in the recent papers (see, 
e.g., Drut & Lahde (2008), Armour et al. (2009)) show that graphene can turn into 
excitonic insulator state while being suspended in vacuum. Another possibility is 
electron-electron pairing in graphene layer, mediated either by phonons, plasmons 
(Uchoa & Castro Neto 2007), or by Coulomb interaction, acting as attractive in cer- 
tain channels (i.e. resonating valence bond mechanisms, proposed by Black-Schaffer 
& Doniach (2007), Honerkamp (2008), or anisotropic electron scattering near Van 
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Hove singularity by Gonzalez (2008)). The other two possible mechanisms of estab- 
lishing a coherent state in graphene are the proximity-induced superconductivity 
(Heersche et al. 2007, Beenakker 2008), and the pairing of spatially separated elec- 
trons and holes in graphene bilayer (Lozovik & Sokolik 2008 o, Min et al. 2008, 
Zhang & Joglekar 2008), analogous to electron- hole pairing in coupled quantum 
wells (Lozovik & Yudson 1975,1976, Lozovik & Poushnov 1997, Lozovik & Bcrman 
1997). 

We consider electron-electron pairing due to optical in-plane phonons, repre- 
sented by two pairs of doubly-degenerate modes (Piscanec et al. 2004, Basko & 
Aleiner 2008), and due to out-of-plane (flexural) acoustical and optical modes; the 
out-of-plane modes interact with electrons quadratically (Mariani & von Oppen 
2008, Khveshchenko 2009). We demonstrate that each phonon mode in graphene 
provides a contribution to effective electron-electron interaction, dependent on both 
its symmetry and the structure of the order parameter with respect to electron val- 
leys. Estimates of the coupling constants show that out-of-plane phonons do not 
cause a pairing with any observable critical temperatures, however the in-plane 
optical phonons can lead to the pairing in heavily doped graphene. 

Electron-hole pairing in graphene bilayer in the weak coupling regime is of BCS 
type, and affects only the conduction band of the electron-doped graphene layer and 
the valence band of the hole-doped layer (Lozovik & Sokolik 2008 a). On increase 
of the coupling strength, the pairing become multi-band, involving also the valence 
band of the electron-doped layer and the conduction band of the hole-doped layer 
(Lozovik & Sokolik 2009, Lozovik& Sokolik 2010). Such "ultrarelativistic" regime 
of pairing occurs due to absence of localized pairs in graphene (Lozovik & Sokolik 
20086, Sabio et al. 2009) — in contrast to usual systems of attracting fermions, 
where crossover to a gas of local pairs at strong coupling occurs (Nozieres & Schmitt- 
Rink 1985). 

The estimates of a critical temperature in graphene bilayer within the framework 
of one-band BCS model with taking into account static screening of electron-hole 
interaction give unobservably small values (Kharitonov & Efetov 2008a,6). However 
the estimates using unscreened interaction, or with the statically screened interac- 
tion, but within the multi-band model, provide much larger values of the critical 
temperature (Min et al. 2008, Zhang & Joglekar 2008, Lozovik & Sokolik 2009). In 
this paper, we consider the strong coupling regime within the framework of multi- 
band model, taking into account a dynamical screening of electron-hole interaction. 
The dynamical effects manifest themselves as virtual plasmons and single-particle 
excitations, which contribution to the interaction is repulsive and thus competes 
with the "direct" Coulomb attraction. 

The article is organized as follows. In Sec. 2 we derive and solve the gap equations 
for the phonon-mediated electron-electron pairing in graphene. In Sec. 3 we study 
electron- hole pairing in graphene bilayer using Eliashberg-type equations. Sec. 3 is 
devoted to conclusions. 

2. Phonon-mediated electron-electron pairing in graphene 

Electrons in graphene populate two interpenetrating triangular lattices A and B, 
composing the bipartite graphene lattice, and two "valleys" K and K' ^ —K 
in momentum space, therefore it is convenient to describe electrons by the ef- 

Article submitted to Royal Society 



Electron- electron and electron-hole pairing in graphene structures 



fective four-component wave function (Castro Neto et al. 2009). Analogously to 
Gusynin et al. (2007), we introduce the four-component electron destruction oper- 
ator ^p = (a^j , -., bj^.p, bj^,,0, a^/ 1 p-)"^, where the operators ap and bp correspond 
to sublattices A and B respectively. In the Heisenberg representation, ^p evolves 
according to the Dirac-type equation: 

Pm7^^p = 0, fi = 0,1,2. (2.1) 

The "covariant" coordinates p° — {i/vp){d/dt), p^'^ = Px,y, Pfj. = {p^,—p^,—p^} 
are used, where up ~ 10^ m/s is the Fermi velocity. The gamma matrices are in the 
Weyl representation: 

The Hamiltonian of linear electron-phonon coupling can be written in the gen- 
eral form (see Fig. H^a)): 

^il-t -^H 4?*.>.^r,* .$,^. (2.2) 

Here g^ and F^ are coupling amplitude and interaction vertex for the /x-th phonon 

mode, ^g^i = Cg^j + c_-. , where Cg'^ is the phonon destruction operator; ^p- = ^^7*^ 
is the Dirac-conjugated spinor and S is the system area. 

We take into account two pairs of degenerate in-plane optical phonon modes, 
most strongly coupled to graphene electrons (Piscanec et al. 2004, Basko & Aleiner 
2008): Al and Bi modes (we denote them by /x = 1, 2) with the momentum q = ±K 
and energy ujj^ sa 0.170 eV, and E2x and E2y modes (/x = 3,4) with g = F, wr ~ 
0.196 eV. The coupling constants (Piscanec et al. 2004) and interaction vertices 
(Basko & Aleiner 2008) for these modes are: gll'^^ « 1.34eV- A, gll'"^ « 0.86eV- A; 

Fi = /, F2 = ^7^ Ta = -7^7', r4 = -7'7'- "' 

The Hamiltonian of quadratic interaction of graphene electrons with out-of- 
plane phonons is more complicated. The major part of the interaction, resulting 
from the deformation potential (Suzuura & Ando 2002, Mariani & von Oppen 
2008), can be written in the form: 



^cnf-|EEEB'«.-+^jv^ 



\_ p p ^ 

ffwffcra' j=l 






V^gi [ 



iqdj pB _ pA 
q{7 qa 



^iq'dj^B _ gA 

q'<T' q'(T' 



AM^/UJ^^UJ^ 



(2.3) 



where g' = p'—p — g, M is the carbon atom mass, 51 ~ 20 — 30eVand dj (j = 1,2,3) 
are the vectors connecting an atom from the A sublatticc with its nearest neighbors. 
The summation over p, p' and q in (j2.3p is performed over the first Brillouin zone 
of graphene; e-J and cu^a are polarizations and frequencies of acoustical {a — 1) 
and optical (a — 2) out-of-plane phonons branches. 
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Figure 1. Diagrammatic representations of electron-phonon coupling Hamiltonians H2.2 
and (|2.4p for linear (a) and quadratic (b) couplings respectively. The corresponding self 
consistent gap equations (|2.6|) and (12.71) are shown in (c) and (d) respectively. 



Rewriting (|2.3p in the four-component spinor notation requires splitting of the 
summation over electron momentums p, p' among the valleys ±K. The result is 
(see Fig. [IJb); see also Lozovik & Sokolik (in press) for details): 
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where q' ~ Q — q. 

Similarly to Pisarski & Rischke (1999), we describe the pairing by the set of ma- 
trix Green functions in Matsubara representation Gij{p,T) ~ — I^T'^^l' (j)'^ s (0)), 

where '^~ = ^I^p-, *l = ^cp = C'^-p' is the charge-conjugated spinor and 
C = ^7^7^ is the charge-conjugation matrix. The anomalous Green functions 6*12 
and G21 are responsible for a Cooper pair condensate. The Gor'kov equations, 
describing the pairing in the mean-field approximation, are: 

G,,{p) = k,Gf{p) + G|°)(p)A,^3-»(p)G3-.,,(p), (2.5) 

where p — {po = iTTT{2k + l),_p}, G[ ^ip) — [l^{po ± A^) — vfJPI^'^ are the free 
particle Green functions, following from (12. ip and /i is the chemical potential in 
graphene. 

The expressions for the anomalous self-energy A21 (the other component A12 = 
7°AJ]^7'') in (12. 5p for the cases of linear (j2.2p and quadratic (12. 4P electron-phonon 
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interaction Hamiltonians are, respectively (see Fig. [ljc,d)), 

T 



A2i(p) = ~^Y.9lD,ip-p')T^G2i{p')r,, (2.6) 



A2i(p) = ^EE^'^(«)^-'(^+p-p'-'?)^S'G2ib')^i^S- (2.7) 



Here the charge-conjugated vertices F^ ~ C "^F^C and V^^^, = C ^V^^J, C are 
introduced, and D^{q) = 2uj^f^/{qQ — w?. ) is the phonon Green function. 

To solve the Gor'kov equations (|2.5p . we assume that the pairing is diagonal 
with respect to the conduction and valence bands, but suppose, that the structure 
of the order parameter with respect to electron valleys is parameterized by arbi- 
trary SU{2) matrix. The band-diagonal order parameter can be represented as a 
decomposition over projection operators 'P±{p) = (1 ± 7'^7p)/2 on conduction and 
valence bands, where p = p/|p| (Pisarski & Rischke 1999). Further, rotation of the 
order parameter in the valley space can be performed by means of three genera- 
tors Ti = 7^, T2 — 7-^7^ and Ta = 17'^, obeying the algebra of the Pauli matrices 
(Gusynin et al. 2007). Thus, the explicit form of A21 is: 

A21 (p) = e-~^ [A+ {p)V+ (p) + A_ {p)V- [p)] , (2.8) 

where A-i-(p) are the gaps in conduction and valence bands, and the three-dimensi- 
onal vector V defines the valley structure of the order parameter. 

Substituting (|2.8p in (j2.6|) and (|2.7I) , we can derive the system of two coupled 
gap equations for A±(p), having the form: 



where E±{p) = \/{vf\p\ T A*)^ + ^±{p) a-i'e the energies of Bogolyubov excitations 

in conduction and valence bands. The effective interband {a = /3) and intraband 
{a — —/3) interactions in the cases of linear and quadratic elcctron-phonon couplings 
are, respectively, 

Kf^iP^P'-^^) = lj29lD,ip-p')Sp [7'„(p)e-»^~r^7"e»^>Mp')7°rJ , (2.10) 



2 



T- 



A'^p{p,p';v) ^ --J2J2^-(i^^-'(Q + P- P' - 1) 

qQ <^<^' 



xSp 



P.(p)e-»-^>£,7°e™~^>,(p')7°V^i,5 



(2.11) 



We can find analytical solutions of (I2.9p - (|2.1ip in the regime of high doping 
of graphene, when ^ is greater than the characteristic phonon frequencies. More- 
over, high doping facilitates the pairing due to larger density of states at the 
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Figure 2. Partial dimensionless Eliashberg functions for out-of-plane phonons. Contribu- 
tions of the two-phonon processes, leaving the electron in its initial valley [Zt, solid line), 
and flipping it into the opposite valley {Zj^, dotted line) are shown. 

Fermi level N = fi/2TTVp. In this case the pairing is effectively one-band, and 
the Eliashberg equations at T = 0, following from (|2.9p . can be derived by setting 
A+(p) = A9(a;o —po), A_ = (similarly to Lozovik et al. 2010): 



1 = 2 



d(jj 



Vi 



w^ 



dv 



to + V 



(2.12) 



Here loq is a cutoff frequency of the order of phonon frequencies, and the Eliash- 
berg functions a'^lv)F{v) for in-plane and out-of-plane phonons respectively can 
be represented as: 



.l{v)F{v) = NalSiv - ujM^\v)+Ugl5{v ^^)R%\v) 



(2.13) 



al{u)F{v) = Zr{u)RP{v) + Z^[u)R't>{v) 



^' 



(2.14) 



The partial Eliashberg functions of out-of-plane phonons can be calculated within 
the simple phonon model (see Lozovik & Sokolik (in press)) and reduced to di- 
mensionless functions: Z^{v) — {8lJ\fg1/2a^APuj'p2)Zfi{x), where x = Qv/ujt2, 
uj-p2 ~ 0.11 eV. The functions Z^{x) are shown in Fig. [51 

Eqs. (|2.13p - (|2.14p show that the Eliashberg functions consist of two parts, corre- 
sponding to phonon-mediated electron-electron interaction processes, leaving both 
electrons in their valleys (F-terms), and flipping both electrons into opposite val- 
leys (iiT-terms). The factors R '- (w), dependent on the valley structure of the order 
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parameter, determine the signs and amplitudes of these contributions: 
i?W ^^>) ^ _ 2Qg2 „ + (1 _ 2wi ) sin^ V, R^V (v) = - cos^ v + vl sin^ v, 

r'^\v) = 1, Rf{v) = {-vl + vj)sm^v. (2.15) 

All the factors (|2.15p except i?p can vary in the range from —1 (the effective 
phonon- mediated repulsion) to +1 (the effective attraction). 

The system under pairing conditions should prefer the valley structure v of the 
order parameter, which provides the maximal gap and thus the maximal a'iiiy) F (ly) . 
For in-plane phonons, the gap, found from (|2.12p - (|2.13p . is A ^ Wp j^exp{ — 1/A}, 
where the effective coupling constant A == 2/^ {dv/v)a'^{y)F{v) consists of the 
partial coupling constants: A = XtRt{v) + Xj^Rg-{v), A^ — 2A/'5^/a;^. When 
A^ > 2Ar, the preferable pairing structure is v = {7r/2, 0, 0}: in this case Rj^ = 1 
(scalar and pseudoscalar phonons give rise to effective electron-electron attrac- 
tion), Rt = —1 (pseudovector phonons cause repulsion). At Xg- < 2Ar, we have 
V = {0, (7r/2)cos(^, (7r/2)sin(^}, and Rr = I (pseudovector phonons cause attrac- 
tion), Rj^ = (contributions from scalar and pseudoscalar modes cancel each 
other). Actually, at high values of dielectric permittivity of surrounding medium, 
A^ < 2Ar; at lower permittivity, the Coulomb interaction renormalizes A^ towards 
higher values, so the relation A^j > 2Ar can be satisfied (Basko & Aleiner 2008). 

For out-of-plane phonons, the preferable valley structure of the order parame- 
ter is V ~ {0,0, 7r/2}, when Rr = Rj^ = 1. Numerical estimates of the coupling 
constant for out-of-plane phonons show very small values (lO"'^ by the order of mag- 
nitude), thus out-of-plane phonons cannot provide any observable electron pairing 
in graphene. However, in-plane optical phonons can provide observable pairing at 
heavy doping of graphene. 

3. Electron-hole pairing in graphene bilayer 

The pairing of spatially separated electrons and holes in graphene bilayer occurs 
due to Coulomb attraction. This longitudinally- vectorial interaction is described by 
the vertex F = 7° within the framework of matrix diagrammatic technique, em- 
ployed in the previous section. Any interaction with such vertex lead to effective 
electron-electron interaction, independent on the valley structure of the order pa- 
rameter. Under band-diagonal pairing, the system of self-consistent equations for 
the conduction- and valence-band gap functions A±(p) is similar to (|2.9p : 

AUp) = -^J2 ^^t^^' ^(\P~^\^Po-Po)mp'): (3-1) 

where Fi3[p') = A^(p')/[pQ^ — E'l{p')\ is the anomalous Green function and V{q,uj) 
is the dynamically screened electron-electron interaction. 

The system (j3.ip can be solved in the spirit of BCS theory (Bardcen et al. 1957), 
i.e. in the static approximation, when one puts a; = in V{q, lo) and assume, that 
^±(^7'^) do not depend on ui and are non-zero in some range of p, corresponding 
to neighborhood of the Fermi surface. Such calculations (Lozovik & Sokolik 2009) 
showed that in multi-band pairing regime the gap depends exponentially on energy 
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width of the pairing region and thus can be very large. Here we go beyond the static 
approximation and take into account the frequency dependence of the screened 
interaction. 

The pairing interaction V{q,uj) can be calculated in the random phase approx- 
imation, well-justified in graphene bilayer due to large number of fermionic flavors, 
equal to 8 (Kharitonov & Efetov 2008a,6; see also Apenko et al. 1982): 

y{(i.^) = ^ _ 2vqIV{q,uj) + «2n2(g, w)(l - e-^iD) ' ^^-^^ 

where Vq = l-ne^ jeq is the bare Coulomb interaction, e is a dielectric permittivity 
of surrounding medium, n((7, w) is a polarization operator of each graphene layer. 
Hereafter we consider the case of small interlayer distance Z?, when ppZ? <C 1, 
Pf = m/'^^f is the Fermi momentum. At w = 0, (13. 2p reduces to the statically 
screened interaction, equal in dimensionless form to MVijq^G) = rsl{qlpF + 8rs), 
where r^ — e- /ev-p sa 2.19/e determines the coupling strength. There exist two 
plasmon branches in the system, corresponding to zeros of denominator of p.2p : 
the lower branch with the dispersion W-(g) ~ v-pq, and the upper branch with the 
square-root dispersion io^{q) « 2/iyVs(g/pF) at small q and almost linear dispersion 
at large q. When w± [q)+v-pq > 2/i, the plasmons in graphene acquire a finite lifetime 
due to interband transitions (Wunsch et al. 2006, Hwang & Das Sarma 2007). At 
cj — > oo, the potential V{q,uj) becomes unscreened: J\fV{q, oo) — Ts/iq/pp)- 

Using the spectral representations oi Fp{p') and V{q) in p.ip and summing over 
Po at T = 0, we get the Eliashberg-type equations: 

A /^ , v^ f dp' 1 + a/3pp' f duj\ ^ ,^, n\ f dv 

p y 

y.luiV{p~p',uj~uj' + i5){— ^ -,-— ^ ^) >. (3.3) 

We assume that A(p, w) is real (it is well-justified near u) = 0), therefore ImF^(p') — 
—TrS{uj' — Efs{p'))Aj3{p')/2Ej3{p'). Then, we assume that the argument of the ^- 
function in this expression vanishes at some unambiguous uj' = Cj{p'). This allows 
us to handle only "on-shell" gap functions and Bogolyubov energies: Aq(|p|) = 
Aa(p, w(p)), Ea{\p\) = Ea{p,oJ{p)). Rewriting (j3.3p in terms of the on-shell quan- 
tities and assuming a = +1, p = pp in its left-hand side, we get: 



A+{pp)^Y. 



dp' l + ppf Apjp') 
(2^)2 2 2Efiip') 



^^-^- + Z f^'^^^P-P''^ + '^^ E,{t') + - ^ ■ ^^-^^ 



TT 

b 



Here we also neglected uj — A+{pp) in the second term in the braces. 

To demonstrate the influence of dynamical screening of V{q, w) on the gap value, 
we will simplify the equation p.4p further. Firstly, we suppose that the both on- 
shell gap functions Ap{p) are equal to each other; this seems plausible at large Tg, 
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Figure 3. On-shell interaction, represented by the expression in the braces of (|3.5|) . in the 
intraband channel (/3 = +1) at rs = 2. Solid line: full on-shell interaction, dash-dotted 
line: unscreened Coulomb interaction only, dotted line: unscreened Coulomb interaction 
with the contribution of undamped plasmons, dashed line: statically screened interaction. 



as discussed by Lozovik & Sokolik 2010, and allows us to neglect the angular factor 
pp' in p.4p . Secondly, we assume A±{p) = A/(p), where A = A+{pp) is the gap 
at the Fermi surface and f{p) is some trial function, equal to 1 at p = pp and 
having the asymptotics f{p) ex 1/p at p — > oo, caused by the leading contribution 
of the unscreened Coulomb interaction in p.3p . For computational purposes, it is 
also convenient to single out the contribution of undamped plasmons of the higher 
branch as: lmy{q,i' + i6) — Imy{q,i') — 7r(5[i^ — uJ-^-{q)]A{q). Here the spectral 
weight of the higher branch plasmons is A{q) = —Vq/2[dll{q, uj) / duj]\i^^^^(^qj when 
q-{-oj+{q) < 2/x and zero otherwise. Thus, the equation p.4p for the gap reduces to: 



1 = 



2^^ 

/3 



dp Aff(p) 
(2^)2 2Epip) 



2A{p- pf) 



^P-PF 



Epip) +uj+{p-pf) 



+ 



OO 

-- / ImV{p-pF,i^) 



dv 



Mp) 



(3.5) 



The expression in braces in p.Sp is the effective on-shell interaction, which 
incorporates the effects of dynamical screening of the pairing interaction. It differs 
from the statically screened interaction V{q,u!), employed by Lozovik & Sokolik 
2009, Lozovik & Sokolik 2010 and is naturally divided into three contributions (see 
Fig. 13]): a) attractive unscreened Coulomb interaction, b) repulsive contribution 
due to virtual undamped higher-branch plasmons, c) repulsive contribution due to 
virtual damped higher-branch plasmons and single-particle intra- and interband 
excitation continuums at a; < vpq and uj + vpq > 2/i respectively (Wunsch et al. 
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Figure 4. Values of the gap at the Fermi surface A, normalized on the chemical poten- 
tial /i, as functions of r-g, calculated with taking into account various contributions to 
the effective on-shell interaction. Solid line: full on-shell interaction, dash-dotted line: un- 
screened Coulomb interaction only, dotted line: unscreened Coulomb interaction with the 
contribution of undamped plasmons, dashed line: statically screened interaction. 



2006, Hwang & Das Sarma 2007). However, strictly on the Fermi surface (at /3 = +1, 
p = Pf) the effective on-shell interaction coincides with the statically-screened one. 

Assuming that the characteristic momentum, at which the on-shell gap A±(p) 
decays, is of the order of the Fermi momentum and employing /(p) = Pf/(|p~Pf| + 
Pf) as a trial function, we can study the influence of dynamical screening of the 
pairing interaction on the gap value A. Fig. S] shows the results of the numerical so- 
lution of p. 51) with taking into account various contributions to the effective on-shell 
interaction. Taking into account only the unscreened Coulomb attraction results in 
huge gap values, reported previously by Min tt al. (2008), Zhang & Joglekar (2008). 
The addition of the undamped plasmons repulsive contribution eliminates the log- 
arithmic singularity of the effective interaction on the Fermi momentum (Fig. |3]), 
but do not change the gap values essentially at large Vg. Finally, the contribution 
of the damped plasmons and and single-particle excitations lowers the gap down to 
the values of the order of 0.01/z at maximum (several Kelvins at /i ^ 0.1 eV). 

When the on-shell interaction is replaced by the statically screened one V{q, 0), 
the gap is unobservably small, if we use the trial function /(p), spread over the 
momentums of the order oi pr] this is in argeement with the BCS-type estimates 
by Kharitonov & Efetov (2008 a, 6). However, as shown by Lozovik & Sokolik (2009), 
the pairing tend to occupy much larger region of the momentums of the order of 
SvsPf and results in large gap values, if we use the statically screened potential 
as the pairing potential. In our case, when we take into account the dynamical 
effects and naturally assume that the on-shell gaps A±{p) decay at p ^ pp, the gap 
turns out to take the values, several orders of magnitude smaller than these with 
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unscreened Coulomb interaction, but, at the same time, several orders of magnitude 
larger than the BCS-type estimates. 



4. Conclusions 

We have considered electron-electron pairing in graphene, mediated by in-plane and 
out-of-plane phonons, and electron-hole pairing in graphene bilayer, mediated by 
the screened Coulomb interaction. In both cases we consider the generally multi- 
band pairing with the s-wave order parameter, diagonal with respect to valence and 
conduction bands of paired particles. Moreover, we take into account the frequency 
dependence of pairing interaction in both cases, deriving and solving two-band 
Eliashberg-type equations. 

The consideration of phonon-mediated pairing is performed by resolution of 
electron-phonon interaction with respect to sublattice and valley degrees of freedom 
of electrons in graphene and taking into account a possibility of different structures 
of the order parameter in valley space. We demonstrate that contributions of dif- 
ferent phonon modes in graphene to effective electron-electron interaction, entering 
the Eliashberg-type equations, depend both on symmetries of these modes and on 
the structure of the order parameter in valley space. 

The coupling of graphene electrons with out-of-plane (flexural) phonon modes is 
quadratic and leads to unusual form of effective electron-electron phonon-mediated 
interaction, which includes integration on frequency and integration on momentum 
over the whole Brillouin zone within the phonon loop. The estimates of effective 
coupling constants show that the pairing due to in-plane phonon modes can occur 
at high doping of graphene, while the pairing due to out-of-plane phonons does not 
occur at observable temperatures. 

Eliashberg-type equations for electron-hole pairing in graphene bilayer, written 
in the present paper in terms of the "on-shell" gap functions, allow to estimate the 
role of dynamical effects. The effective on-shell dynamically screened interaction, 
entering the gap equations, can be represented as a sum of attractive unscreened 
Coulomb interaction, repulsive contribution due to undamped virtual plasmons and 
combined repulsive contribution of the damped plasmons and of the continuum of 
single-particle excitations. 

The unscreened Coulomb interaction on its own provides large values of the 
gap, which are only slightly reduced with taking into account the undamped plas- 
mons. Inclusion of the damped plasmons and single-particle excitations lower down 
the gap by several orders of magnitude. This result demonstrates the significant 
competition between the bare Coulomb attraction and virtual excitations in the 
system, responsible for its dynamical screening. However, the estimates of the gap, 
calculated with taking into account the full on-shell potential, are by several orders 
of magnitude larger than the BCS-type estimates and can reach several Kelvins at 
strong coupling. 
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